Pore scale mixing and macroscopic solute dispersion regimes in polymer flows inside 

2D model networks. 
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A change of solute dispersion regime with the flow velocity has been studied both at the macro- 
scopic and pore scales in a transparent array of capillary channels using an optical technique allowing 
for simultaneous local and global concentration mappings. Two solutions of different polymer con- 
centrations (500 and 1000 ppm) have been used at different Peclet numbers. At the macroscopic 
scale, the displacement front displays a diffusive spreading: for Pe < 10, the dispersivity Id is 
constant with Pe and increases with the polymer concentration; for Pe > 10, Id increases as Pe^'"^^ 
and is similar for the two concentrations. At the local scale, a time lag between the saturations of 
channels parallel and perpendicular to the mean flow has been observed and studied as a function 
of the flow rate. These local measurements suggest that the change of dispersion regime is related 
to variations of the degree of mixing at the junctions. For Pe < 10, complete mixing leads to 
pure geometrical dispersion enhanced for shear thinning fluids; for Pe > 10 weaker mixing results 
in higher correlation lengths along flow paths parallel to the mean flow and in a combination of 
geometrical and Taylor dispersion. 

PACS numbers: 47.55.Mh, 7.25. Jn, 05.60+w 



I. INTRODUCTION 

The problem of solute transport in porous media is 
relevant to many environmental, water supply and in- 
dustrial processes [jj, l2|. In addition, tracer dispersion is 
a useful tool to analyse porous media heterogeneities ^ . 
Solute dispersion is often measured by monitoring solute 
concentration variations at the outlet of the sample fol- 
lowing a pulse or step-like injection at the inlet: Under- 
standing fully the dispersion mechanisms requires how- 
ever informations on the local concentration distribution 
and on local mixing at the pore scale. These goals can be 
reached by using NMR-Imaging, CAT-Scan or acoustical 
techniques [3| : however, such techniques are either costly 
or have a limited resolution. In addition, they often put 
strong constraints on the characteristics of the fluid pairs. 
We have chosen instead in the present work to use a 2D 
transparent square network of ducts of random widths 
allowing for easy visualizations of mixing and transport 
processes at the local scale. Similar systems were previ- 
ously used successfully to investigate two phase flows [1] 
and miscible displacements of Newtonian fluids . 
The key feature of the present experiments is to com- 
bine macroscopic and local scale measurements to esti- 
mate the influence of pore scale processes on dispersion. 
For this purpose, dye is used as a solute and high res- 
olution maps of the relative concentration distribution 



are obtained through calibrated light absorption mea- 
surements. This allows one to determine quantitatively 
the global dispersion coefficient D and the dispersivity 
Id while measuring simultaneously the time lag between 
the invasion of channels parallel and transverse to the 
mean flow. Further informations are also obtained from 
the complex geometry of an isoconcentration line. 
This approach have allowed us to observe a change of 
dispersion regime for a Peclet number of the order of 10 
while the local measurements suggest interpretations in 
terms of variations of the degree of mixing at the junc- 
tions. This confirms previous suggestions 0, on the 
influence of mixing at the pore scale on macroscopic dis- 
persion. 

Another feature of our experiments is the use of shear 
thinning polymer solutions of the type encountered in 
many industrial processes in petroleum, chemical and 
civil engineering [9| . In addition to these practical appli- 
cations, previous measurements at the laboratory scale 
on glass bead packings have shown [l^, Uli a significant 
enhancement of tracer dispersion compared to the New- 
tonian case. This enhancement depends on the polymer 
concentration and will represent a useful additional input 
for our interpretations. 



II. DISPERSION MECHANISMS IN 3D 
POROUS MEDIA AND 2D NETWORKS 
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In homogeneous SD porous media, the macroscopic 
concentration C of a tracer (i.e. averaged over a rep- 
resentative elementary volume) satisfies the convection- 
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diffusion equation : 

dC{x,t) _ dC{x,t) d^C{x,t) 

dt dx ^ dx^ ^' 

where D is the longitudinal dispersion coefficient, U the 
mean velocity of the fluid (parallel to x) and C is as- 
sumed to be constant in a section of the sample normal 
to U. In the following, equation ^ is shown to be also 
valid for the transparent model used in the present work. 
The value of D is determined by two main physical mech- 
anisms : molecular diffusion and advection by the com- 
plex velocity field inside the medium (the local flow veloc- 
ity varies both inside individual flow channels and from 
one channel to another). The relative order of magni- 
tude of these two effects is characterized by the Peclet 
number: Pe = Ua/2 Dm {Dm is the molecular diffusion 
coefficient and a a characteristic length of the medium, 
here the average channel width). 

Various dispersion regimes are observed in usual porous 
media [1,]. At very low Peclet numbers (Pe < 1), molec- 
ular diffusion is dominant and smoothes out local con- 
centration variations. 

At higher Pe values, the distribution of the channel 
widths induces short range variations of the magnitude 
and the direction of the local velocity. One can then 
consider that tracer particles experience a random walk 
inside the pore volume with a velocity varying both in 
magnitude and direction relative to the mean flow veloc- 
ity U. The typical length / of the channels represents 
the length of the steps and their characteristic duration 
T is T « l/U . A classical feature of random walks is that 
the corresponding diffusion coefficient (here equal to D) 
satisfies D cx P/r = Ul. The proportionality constant 
depends both on the disorder of the medium and on the 
rheology of the fluid. In this so-called geometrical disper- 
sion, the coefficient D should then be proportional to U . 
The 2D networks of interest in the present work have 
several specific features and D should be very much in- 
fiuenced by the redistribution of the incoming tracer be- 
tween the channels leaving a junction. This redistribu- 
tion strongly depends on the local structure of the flow 
field and on the Peclet number [l^, [l^l ■ At low Pe val- 
ues, the transit time through a junction is large enough 
for tracer to cross streamlines by molecular diffusion and 
one may assume a perfect mixing. In the other limit 
Pe » 1, molecular diffusion is negligible: the path of 
the tracer particles coincides with the flow lines and is 
determined by the flow field in the junctions and by the 
location of the particles in the fiow section. 
In the case of 2D networks with small variations of the 
channel apertures, the flow field is close to that in a pe- 
riodic square channel with a mean fiow parallel to one 
of the axis: the major part of the fiow is localized in 
longitudinal channels parallel to the axis where the ve- 
locity is high while that in transverse channels is small. 
Tracer particles remain then inside sequences of chan- 
nels parallel to the mean fiow for a long distance without 
moving sideways. Taylor-like dispersion p^ . [l5| similar 



to that encountered in capillary tubes, between parallel 
plates or inside periodic structures may then develop: it 
corresponds to a balance between (a) spreading due to 
velocity gradients between the centers of the fiow chan- 
nels and their walls and (b) molecular diffusion across 
the flow lines. The increased dispersion in 2D periodic 
networks when flow is parallel to one of the axis [t'I may 
for instance result from a combination of geometrical and 
Taylor dispersion. 

Another issue of the present work is the influence of the 
shear-thinning properties of the fluids. They may influ- 
ence the dispersion process in different ways: on the one 
hand, when the viscosity /i decreases with the shear rate 
7, the velocity profile in individual channels becomes fiat 
in the center of the channel. In simple geometries like 
capillary tubes, the corresponding Taylor dispersion co- 
efficient D is then lower than for a Newtonian fluid al- 
though one has still D oc Pe^ [T6l|. On the other hand, 
numerical investigations suggest that the flow of shear 
thinning fluids is localized in a smaller number of flow 
paths than for Newtonian fluids [13, As a result, 

geometrical dispersion reflecting the distribution of the 
local velocities should increase as is indeed observed ex- 
perimentally [13, EH ■ A key point is here the influence of 
the fluid rheology on tracer mixing at the intersections 
between channels: to our knowledge, no previous experi- 
mental or numerical work has dealt with this issue which 
is an important point. 



III. EXPERIMENTAL SET-UP AND 
PROCEDURE 



A. Experimental models and injection set-up 

The model porous medium is a two dimensional square 
network of channels of random aperture: these models 
have been realized by casting a transparent resin on a 
photographically etched mold as described in ref . [^ . The 
model contains a square network of 140 x 140 channels 
with an individual length equal to I = 0.67 mm and a 
depth of 0.5 mm; the width follows a discrete, log-normal 
distribution with 7 values between 0.1 and 0.6 mm (the 
average width is a = 0.33 and its standard deviation 
cr(a) = 0.11 mm). The mesh size of the network is 
equal to 1 mm. The overall size of the model is 150 x 
140 mm and the two facing lateral sides are sealed. On 
the two others, the channels are directly connected to the 
outside. The total pore volume is close to 6.09 x 10'^ mm^. 
Following the definition of Bruderer et al [l^ , the degree 
of heterogeneity of the network can be caracterized by 
the normalized standard deviation a{a)/a. In the present 
work : a{a)/a ~ 0.33. 
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FIG. 1: Experimental setup for miscible displacement mea- 
surements in 2D micromodels. 



that rjoo ~ 10~^ Pa.s (the value for the solvent i.e. wa- 
ter). In Eq.(I21), 7o corresponds to a crossover between 
two regimes. On the one hand, 7 < 70, the viscosity 
r\ tends towards 770, and the fluid displays a Newtonian 
behavior ("Newtonian plateau" regime). On the other 
hand, for 7 S> T'o, the effective viscosity decreases with 
the shear rate following the power law oc 'y("~i) (n = 1 
for a Newtonian fluid.) It should finally be noted that 
the high effective viscosity of these solutions at low shear 
rates avoids the appearance of buoyancy induced insta- 
bilities at low shear rates and helps stabilize the fluid 
displacement. 



Polymer Cone, 
ppm 

500 
1000 



70 



mPa.s 



0.38 ± 0.04 0.077 ± 0.018 410 ± 33 
0.26 ± 0.02 0.026 ± 0.004 4500 ± 340 



B. Fluid preparation and characterization 



TABLE I: Rheological parameters of scleroglucan solutions 
used in the flow experiments. 



The fluids used in the experiments are shear thinning 
water-polymer solutions. The shear thinning fluids are 
solutions of 500 to 1000 ppm of high molecular weight 
scleroglucan in high purity water (Millipore - Milli-Q 
grade). Scleroglucan (Sanofi Bioindustries, France) is a 
polysaccharide with a semi rigid molecule (persistence 
length ~ 180 x 10"^ m); it has been selected because 
it is electrically neutral and its characteristics are there- 
fore independent of the ion (and dye) concentration. All 
solutions are protected from bacterial contamination by 
adding 0.2 g/1 of NaN^. In all experiments, the in- 
jected and displaced fluids are identical except for a small 
amount of Water Blue dye which is added to one of the 
solutions at a concentration of 200 ppm by weight to al- 
low for optical concentration measurements. 
The molecular diffusion coefficient of the dye is deter- 
mined independently in pure water by means of Tay- 
lor dispersion measurements performed separately in a 
capillary tube (the measured value of Dm was close to 
6.5 X lO"** mm^.s"^.) 

The rheological properties of the scleroglucan solutions 
have been characterized using a Contraves LS30 Couette 
rheometer with shear rates 7 ranging from 0.016 s~^ up 
to 87 s^^. The rheological properties of the solutions have 
been verifled to be constant with time within experimen- 
tal error (over a time lapse of 3 days) and identical for 
dyed and transparent solutions with the same polymer 
concentrations. The variation of the effective viscosity 77 
with 7 was then adjusted by a Carreau function : 
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-(?70 - ??oo) + ??o 



(2) 



The values of these different rheological parameters for 
the polymer solutions used in the present work are listed 
in Table HI determining r]^ would have required values 
of 7 outside the measurement range so that we assume 



C. Fluid injection and flow visualization 

The model is placed vertically with its open sides hori- 
zontal : the upper side is fitted with a leak tight adapter 
allowing one to suck fiuid upwards. The lower open side 
is initially slightly dipped into a bath of one of the liq- 
uids and is saturated with this fluid by pumping it slowly 
upwards. After switching off the pump, the bath can 
be lowered until it does not touch any more the model 
(Fig. [T]). The bath is then completely emptied, refilled 
with the other fluid and raised to its initial position. Fi- 
nally, the first fluid is sucked upwards at the upper end 
of the model by the syringe pump. This allows one to 
obtain a front of the displacing fluid which is initially 
perfectly straight. The lower bath rests upon computer 
controlled electronic scales for monitoring the amount of 
fluid which has entered the model. The flow rates used 
in the experiments correspond to mean front velocities 
between 0.005 and 2.5 mm.s^^. 

The model is illuminated from the back by an electronic 
light panel and images are acquired by a 12 bits, high 
stability, digital camera with a 1030 x 1300 pixels reso- 
lution and then recorded by a computer. The pixel size 
is 0.14 mm or 0.4 times the mean width of the chan- 
nels: this allows one to discriminate between the vari- 
ous regions of the pore space. Typically 100 images are 
recorded for each experimental run at time intervals from 
2.5 to 700 s. 



D. Image analysis procedure 

The images are then translated into maps of the rela- 
tive concentration using the following procedure. First, 
a calibration curve is obtained from images of the model 
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saturated with 7 solutions of increasing dye concentra- 
tion C starting from zero up to the concentration used 
in the experiments. The logarithm Ln{I{C) / lo) of the 
transmissivity is then plotted as a function of C; I{C) 
and lo are averages of the light intensity over the model 
and correspond respectively to dye concentrations equal 
to C and 0. Due to the non linear absorbance effect [20|, 
a better fit is obtained with a third order polynomial 
variation of Ln{I{C) / lo) with C than with the linear 
dependence corresponding to Lambert's law. This cali- 
bration is performed every time the locations of the light 
source and of the models are significantly changed. 
For all experiments, reference images are recorded both 
with the micromodel initially saturated with the dis- 
placed fluid and, at long times, when it is fully saturated 
with the injected one. After the fluid displacement has 
been performed, the local concentrations are determined 
pixel by pixel for each image by means of the calibra- 
tion curve. This operation is performed only on pixels 
belonging to channels; non flowing domains are not con- 
sidered. Finally, maps of the local relative concentration 
of the two fluids are obtained by normalizing the local 
concentration between its values in the initial and final 
images. 




FIG. 2: Relative concentration maps for experiments us- 
ing water-scleroglucan solutions of respective concentrations 
500 ppm (a,b) and 1000 ppm (c,d) for flow rate values 
Q = 0.075 ml/mn (b,d) and Q — 1.5 ml/mn (a,c). Grey 
levels code used in the figure : darkest shade — pure dis- 
placed fluid, intermediate - pure displacing fluid - lightest — 
intermediate concentrations 



IV. EXPERIMENTAL RESULTS 



B. Quantitative dispersion measurement procedure 



Qualitative observations of miscible 
displacements 



Figure [2] compares displacement experiments realized 
at a same flow velocity for the two water-scleroglucan 
solutions. For a given solution, narrow structures of the 
front with a lateral extent of a few channel widths appear 
when the velocity increases and reflect local high or low 
velocity zones. For a same flow rate, images obtained 
with the two solutions are qualitatively similar: the size 
parallel to the flow of medium scale front structures (with 
a width of the order of one tenth that of the model) is 
however larger for the more concentrated solution both 
at low and high velocities. 

The stability of the displacement with respect to buoy- 
ancy driven instabilities has also been verified by com- 
paring experiments using the same pair of fluids and ex- 
changing the injected and the displaced fluid jTlj: no 
quantitative difference was measured between the two 
configuraton and no fingering instabilities appeared in 
the unstable configuration. 

While fluid displacement images provide informations on 
the front geometry down to a fraction of the channel size, 
we shall see next that they also allow to determine macro- 
scopic parameters characterizing the process. 



Quantitatively, the global displacement process is ana- 
lyzed from the variations with time and distance parallel 
to the flow of a macroscopic concentration C{x,t): C 
is the average over an interval Ay perpendicular to the 
flow of the local relative concentration C for individual 
pixels located inside the pore volume. The width Ay is 
large enough to average out local fluctuations and small 
enough avoid the influence of the side walls: C was found 
to be independent of Ay within experimental error when 
these conditions are verified. Also, it will be seen below 
that the results are hardly different when Ay corresponds 
to only one channel. 

Figure [3] displays a typical variation of C with time. 
For step-like initial concentration variations correspond- 
ing to our experiments, the solution of the convection- 
dispersion equation ^ is: 



C 



1 - erf 



t - t 



(3) 



As seen in Figure [3l a very good fit of the experimental 
data with this solution (continuous line) is obtained by 
adjusting the two parameters of the equation, namely the 
mean transit time t and the ratio D/U^ (= At^ /2t where 
At^ is the centered second moment) 
Figure [?] displays the variation of the mean transit time 
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t[s] 400 



FIG. 3: Normalized concentration variation C{x,t) as a func- 
tion of time for a displacement experiment using a 500 ppm 
polymer solution with Q = 3.75 ml/mn. C: averaging inter- 
val Ay ~ 35 mm located in the central part of the model. 
Distance from inlet: a; = 80 mm. Continuous line : fit by 
Eq. ®. 
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X [mm] 

FIG. 5: Variation of D with the distance X for dyed fluid 
displacing clear fluid (Polymer concentration : 500 ppm, Q — 
3.75 ml/min, Pe — 317.) - Inset : a close-up at the scale of 4 
channels. 
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curves have been obtained for dyed fluid displacing clear 
fluid, comparison experiments have been realized with 
clear fluid displacingg dyed fluid; no systematic differ- 
ence between the two sets of data was observed, confirm- 
ing that there are no buoyancy driven instabilities. 
The oscillations of t and D as X varies are closely re- 
lated to mixing at the junctions and to the exchange of 
tracers between the transverse channels and the rest of 
the flow (a quantitative analyzis will be presented in sec- 
tion |TVDl) 

The above analysis has been performed for all experi- 
ments realized with both polymer concentrations. In the 
following section, variations with Pe of the dispersion 
coefficient D measured in this way are discussed 



FIG. 4: Variation of t with the distance X (polymer concen- 
tration : 500ppm, Q = 3.75ml /min, Pe = 317) for dyed fluid 
displacing clear fluid. Dashed line : clear fluid displacing dyed 
fluid. Inset : close-up at the scale of 4 channels. Dashed line: 
linear regression of the variation of t with X. 



t with the distance X from the inlet: t increases overall 
linearly with distance, indicating that the mean velocity 
U is constant and can be determined by a linear regres- 
sion of the data. The inset (magnified view of the curve) 
shows that t{X) oscillates about the mean trend (dotted 
line) corresponding to the linear regression. These pe- 
riodic variations are directly related to the structure of 
the network and will be discussed in section HVDI Some 
similar features are observed on the variation of the dis- 
persion coefficient D plotted in figure [5] as a function of 
the distance X from the injection [D is computed from 
the value of D/U^ given by the fit with U equal to X/t.) 
This time, D{X) is globally constant with X; like t, it 
displays periodic oscillations related to the structure of 
the network which will be discussed below. While these 



C. Flow velocity dependence of dispersion 
coefficient 

The variations with the Peclet number Pe of the dis- 
persion coefficient D determined as explained in the pre- 
vious section are displayed in the inset of figure [6] For 
both polymer concentrations, D always increases with 
Pe but two different variation regimes are visible. 
For Pe > 10, values of D corresponding to the two poly- 
mer concentrations fall on top of each other and increase 
like D ^ Pe^ with (3 ~ 1.35. This is in good agreement 
with numerical simulations [l^ realized for a similar ge- 
ometry and degree of heterogeneity (as characterized by 
a{a)/a = 0.33) and for which a power law variation 
with an exponent of the order of 1.3 is also obtained. The 
value significantly higher than 1 of this exponent cannot 
be accounted for solely by geometrical dispersion (which 
would give a value of 1) even if a logarithmic correction 
factor such that D cx Pe LogPe is introduced [2l|. A 
likely hypothesis is that this variation reflects a crossover 
from geometrical dispersion {D cx Pe) to Taylor disper- 
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1000 



FIG. 6; Variation of the dispersivity la (mm) with the 
Peclet number for experiments with water-polymer solutions 
: (□),(■) : 500 ppm concentration - (o), (•) 1000 ppm. Open 
(resp. grey) symbols : averaging interval : 35 (resp. 0.4) 
mesh sizes - Inset : variation of normalized dispersion coeffi- 
cient D / Dm as function of the Peclet number for two polymer 
concentration (same symbols as in main graph). Solid line : 
power law fit for Pe > 10 (exponent 1.35 ± 0.03). 



sion [D oc Pe^): in that range of Peclet numbers, both 
mechanisms would then contribute to the dispersion pro- 
cess. 

For Pe < 10, on the contrary, values of D obtained for 
the 1000 ppm polymer solution are higher than those 
obtained with the 500 ppm one. The variation of D 
with Pe in that range has been studied more precisely 
by plotting the dispersivity Id = D/U in the main graph 
of figure [HI The value of la is approximately constant for 
Pe < 10 for both polymer concentrations and is larger 
for the 1000 pmm solution [l^ = 1.71 mm) than for the 
500 ppm one (Id = 0.8 mm). This implies that, in this 
range of Peclet numbers, the geometrical mechanism con- 
trols dispersion and that the corresponding dispersivity 
increases with the shear thinning character of the fluid 
(at higher Pe values. Id is, in contrast, almost identical 
for the two solutions). 

The values plotted in figure have been obtained by 
averaging the concentration C over 35 channels in the 
direction perpendicular to the mean flow. In order to 
estimate the influence of the heterogeneities of the net- 
work, we performed an analysis in which the concentra- 
tion C is only averaged over Ay = 4 pixels (or about 
0.55 mesh sizes of the lattice): the coordinates y of these 
measurement lines are chosen so that all pixels are inside 
longitudinal channels or junctions and the dispersivity is 
determined as above. 

These values of Id are represented as grey symbols in fig- 
ure [S] and are only slightly lower than those obtained for 
Ay — 35 mm. This shows that there are no large scale 
heterogeneities of the network, such as high or low per- 
meability channels of width significantly larger than the 
mesh size which would increase the dispersivity. 
Since all data points correspond to Pe > 1, the direct in- 



fluence on longitudinal dispersion of molecular diffusion 
is negligible: It has however a strong indirect influence 
at the lower Peclet numbers investigated (Pe < 10). The 
transit time of the tracer inside the junctions or individ- 
ual flow channels is then large enough so that molecular 
diffusion across the flow lines is signiflcant : this influ- 
ences strongly the redistribution of the incoming tracer 
between channels leaving each junction 12]. In the next 
section, we show that, in addition to the determination 
of macroscopic parameters like U , D (or Id), the concen- 
tration maps allow to investigate mixing processes at the 
pore scale or even below. 



D. Tracer exchange dynamics between transverse 
and longitudinal channels. 

The variations with the distance X from the inlet of 
both the mean transit time t{X) (Fig. 2]) and the dis- 
persion coefficient D{X) (Fig. [5]) display periodic oscil- 
lations about respectively an increasing linear trend and 
a constant value. In the following, the oscillations of D 
will be characterized quantitatively by the difference be- 
tween D{X) and its mean value; the variations of t{X) 
will be characterized by its difference with the linear re- 
gression ligne over all data points reflecting the mean 
front velocity U. At a given distance X the time cor- 
responding to this regression is equal to X/U. The dif- 
ference t{X) — X/U is negative when the line located at 
the distance X from the inlet, and over which C{X,t) 
is computed, contains only longitudinal channels paral- 
lel to the flow; it is positive when the line contains both 
transverse channels and junctions. Regarding D[X), it 
is also lower than the mean value when the line X — est. 
contains only longitudinal channels and higher when it 
contains both longitudinal and transverse channels. The 
variations of t{X) reflect the different influence of longi- 
tudinal and transverse channels on transport. As already 
discussed in section|lIl in a weakly disordered square net- 
work like the present one, most convective flux is local- 
ized inside the longitudinal channels (parallel to the mean 
flow) . The mean velocity inside them is then significantly 
higher than in the transverse channels and they get satu- 
rated faster with the displacing fluid. There is therefore 
a time lag between the saturation of the transverse and 
longitudinal channels at a same distance X from the inlet 
which explains the oscillations oit{X) in Figure [H 
Moreover, the respective amplitudes of the successive 
minimas and maximas of ^.{X) — X/U are found to be 
almost constant from one to the other. In the following, 
the time lags will therefore be characterized by the re- 
specting averages 5ti and 5tjt over all minimal and max- 
imal values of t{X) — X/U. 

In the limit of a perfect mixing at the junctions (Pe < 
10), fluid particles do not retain the memory of their 
past trajectory (ie whether they got previously trapped 
inside slow transverse channels). Therefore the time lag 
for X = est lines containing transverse channels should 



7 



reflect directly the residence time in an individual (slow) 
transverse channel, weighted by the volume fraction cor- 
responding to these channels. As long as Pe ^ 1 and 
molecular diffusion is negligible, the residence time in a 
given channel will be inversely proportional to the local 
velocity; the latter is, in turn, proportional to the mean 
velocity U as long as (as in the present case) the Reynolds 
number is low enough and the linear Stokes equation is 
approximately applicable. As a result, the local velocity 
in a channel is proportional to the mean velocity U and 
the time lag Stjt > should vary as 1/U (for similar rea- 
sons Sti < should also vary as l/U). Figure [7] displays 
the variations of Stjt and SU with 1/U for the two poly- 
mer solutions investigated : in both cases, the variation 
is indeed linear for 1/U > 40 s.mm~^ (corresponding to 
Pe < 6). 

For 1/U < 40 {Pe > 6), Stjt (resp. Sti) is higher 
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ties of the fluid : the effective viscosity of the solutions 
increases much more with the polymer concentration in 
slow transverse channels (where the shear rate is low) 
than in fast longitudinal ones : as a result, the contrast 
between the velocities (and therefore the residence times) 
in the longitudinal and transverse channels is enhanced, 
leading to the observed increase of \Stjt\ and \Sti\. This 
enhancement of velocity contrasts for shear thinning flu- 
ids is discussed in more detail in section IVl 
At the opposite limit of low velocities such that Pe < 1 
{1/U > 250 sxmm~^), longitudinal diffusive transfer be- 
comes significant. The increase of the residence times 
with 1/C/ is then limited by molecular diffusion to a value 
of the order of a few P /Dm ('^ 1000s): the variations of 
Stjt and Sti should then level off at high 1/U values. The 
lowest values of Pe are however still too high in our ex- 
periments to observe this effect. 

These results suggest therefore that the type of local ob- 
servations reported here provide important informations 
on mixing processes at the pore scale. As pointed out 
recently [3, [1], these processes may, in turn, influence 
strongly mass transfer at the macroscopic scale. In a sim- 
ilar perspective, we shall now investigate the dependence 
of the geometry of the iso concentration fronts on the 
flow velocity and the polymer concentration. We shall 
also discuss their relation to local mixing in the pores. 



E. Geometry of tracer displacement fronts 



1/Pe 



FIG. 7: Variation of Sti (filled symbols) and Stjt (open sym- 
bols) with 1/Pe. (□,■): experiments realized with 500 ppm 
polymer solutions, (o, •): 1000 ppm. Dashed lines are guides 
for the eyes. 



(resp. lower) than the values extrapolated from the lin- 
ear trend for 1/U > 40 : the transition is observed at the 
same mean velocity for the two solutions. This increase 
of the absolute values of the time lags reflects very likely 
the breakdown of the assumption of perfect mixing at 
junctions or inside individual channels : at high Peclet 
numbers, a solute particle may indeed flow through sev- 
eral junctions and channels without moving across flow 
lines through transverse molecular diffusion. If the lat- 
tice is weakly disordered, the solute may then remain for 
a longer time inside a sequence of longitudinal high veloc- 
ity channels than in the case of a perfect mixing at the 
junctions : the corresponding value of Sti will then be 
lower. Similarly a solute particle may remained trapped 
for a longer time, in slow zones than if mixing was more 
effective at the junctions so that Stjt is increased. 
Although the transition between the above two regimes 
takes place at the same Peclet number for both solutions, 
the absolute amplitude of the variations of Stjt and Sti 
with 1/U is signiflcantly larger for the 1000 ppm one. 
This is direct consequence of the shear thinning proper- 



FIG. 8: Isoconcentration fronts at a mean velocity U = 
0.005 mm/s (Pe — 1.3) for polymer solutions of concentra- 
tions 1000 ppm (upper curve) and 500 ppm (lower curve). 
Mean distance of tracer from inlet : 0.5 x L {L = model 
length). Flow is upward on the figure. Front widths : 
a = 4.5 mm (1000 ppm) and a — 2.6 mm (500 ppm). 

In the present experiments, the pixel size in the concen- 
tration maps is 0.4 times the mean channel width. This 
allows for a study of the tracer distribution in the mixing 
zone at length scales varying from the size of the network 
down to a fraction of the pore size. For practical reasons, 
we shall not use the full spatial concentration distribution 
in the following analysis : we chose instead to character- 
ize its spatial heterogeneity by the isoconcentration lines 
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FIG. 9: Isoconcentration fronts for a mean velocity U = 
0.125 mm/s {Pe = 32) for polymer solutions of concentrations 
1000 ppm (upper curve) and 500 ppm (lower curve). Mean 
distance of tracer from inlet : 0.5 x L {L — model length). 
Flow is upward on the figure. Front widths a — 6.1 mm 
(1000 ppm) and a — 4.6 mm (500 ppm). 



c = 0.5, which are assumed to reflect the displacement 
front geometry: examples of such fronts determined by a 
thresholding procedure are displayed in figures [8] and [9] 
for Pcclet numbers Pe respectively equal to 1.3 and 32. 
The vifidth of the front parallel to the mean flow is larger 
for the more concentrated solution and it increases with 
the Peclet number. Also, at high Pe values, large spikes 
are visible while the front is relatively smooth at lower 
ones. In spite of these differences, the main geomet- 
rical features of the front are similar: large peaks and 
troughs are generally located at the same points for dif- 
ferent flow velocities and polymer concentrations. This 
confirms that irregularities of the front structure are as- 
sociated to deterministic features of the velocity field and 
not to uncontrolled imperfections of the injection. 
Quantitatively, the effective width of the front parallel 
to the flow is characterized in the following by the stan- 
dard deviation a{x) of the distance x of its points from 
the inlet: a{x) satisfies a{x) —< {x — x{y)Y >^^^ hi 
which X is the mean value of x. In figures [5] and [51 x is 
equal to half the length of the model and the values of 
a corresponding to the curves displayed are listed in the 
captions. Figure [TUl displays the variation of cr as a func- 
tion of the mean flow velocity U for x = L/2. For both 
solutions, the width a increases logarithmically with U: 
the value of a is larger for the 1000 ppm solution while 
the slope of its variation with U in figure [TO] is slightly 
lower. 

The values of a are different for the two solutions be- 
cause the effective viscosity decreases faster with the ve- 
locity gradients for the 1000 ppm solution than for the 
500 ppm one. The ratio between the effective viscosities, 
and therefore the velocities in the longitudinal and trans- 




FIG. 10: Effective front width cr as a function of the mean 
front velocity U (mm/s) for a mean distance x = 0.5 x L and 
polymer concentrations equal to 1000 ppm (o) and 500 ppm 
(□). Solid (resp. dotted) lines : fit of the experimental data 
for a 1000 ppm (resp. 500 ppm) polymer solution by equation 
a = a + blog{U) with a = 6.9 ± 0.1 and 6 = 0.4 ± 0.05 (resp. 
a = 6.2 ± 0.2 and 6 = 0.5 ± 0.1. 



verse channels is therefore higher and the front width 
which is directly related to this ratio increases. 
The second major feature of the above results at high 
velocities is the large amplitude of the peaks and troughs 
observed on the front while they are smaller and narrower 
at lower velocities. This, too, may be explained by the re- 
duced tracer mixing at junctions at high Peclet numbers 
(section [IV Dl) : solute remaining for a long distance in- 
side longitudinal, high velocity, channels will contribute 
to the spikes while that moving through a sequence of 
slow lateral channels will contribute to the troughs. At 
lower Peclet numbers, mixing is more efficient and solute 
particules sample more effectively the velocity distribu- 
tion: this reduces the dispersion of the transit times and, 
therefore the amplitude of the peaks and troughs. 



V. DISCUSSION AND CONCLUSIONS 

The local analysis of the transit times and of the front 
geometry have therefore provided important informa- 
tions on mixing inside junctions and flow channels and on 
its dependence on the Peclet number. This information 
greatly helps one to interprete the macroscopic dispersiv- 
ity measurements of section IIVBI Some of the features 
observed are specific to 2D systems while others can oc- 
cur in usual 3D media. 

At low Peclet numbers (typically Pe < 10), the dis- 
persivity Id remains constant with Pe and is lower for 
the 500 ppm polymer solution than for the 1000 ppm 
one. In sections IIV Dl and IIVEI we have seen that, in 
this regime, transverse mixing in junctions and channels 
is very effective so that the correlation length of the mo- 
tion of solute particles is of the order of the length / of 
individual channels. As a result, this motion may be 
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described as a sequence of random steps of varying dura- 
tions and directions inside the medium; this is the geo- 
metrical dispersion regime discussed in section and for 
which D oc Pe (Id = est.). 

In this regime, the factor of two difference of the dis- 
persivities for 500 and 1000 ppm solutions is likely due 
to enhanced velocity contrasts between the fast and slow 
flow regions. It is known, for instance |17j, that the mean 
flow velocity inside a cylindrical channel under a given 
pressure gradient varies as the square of the radius a for 
a Newtonian fluid and as a^+" for a shear thinning fluid 
verifying equation Let us assume that the pressure 
gradient between the ends of flow channels in parallel is 
constant. If (T(a) <C a, the standard deviation 5U of the 
mean velocities in the different channels should scale like 

SU I + n a(a) 

-J7- (4) 

U n a 

where a {a) /a is the normalized standard deviation of the 
channel aperture (see section UlI A[) . Still using the same 
simplistic approach, the typical standard deviation St of 
the transit time along a channel of length / = 0.67 mm 
should be 1/5U. Estimating the dispersion coefficient D 
from the relation D P /6t provides the order of magni- 
tude of the dispersivity : 



n a 

Since n decreases with the polymer concentration. Id 
should therefore increase for a fixed aperture fluctuation 
a{a)/a. Using in equation (O the values of Z, a{a)/a 
and n corresponding to the present experiments (table |l| 
leads to Id — 0.7 mm and 1 mm respectively for the 
500 ppm and 1000 ppm solutions. These estimations 
are close to the experimental values Id = 0.8 mm and 
1.7 mm reported in section HV CI for the same solutions 
(figure [S]). The difference may be due to the assumption 
of identical pressure gradients on different parallel chan- 
nels used to obtain equations 2] and [51 
At higher Peclet numbers Pe > 10, Id is no longer con- 
stant but increases with Pe. This reflects the transition 
towards a second dispersion regime in which mixing is less 
effective. One must then take into account the stretching 
of dye parallel to the flow by local velocity gradients in 
the flow section (dye moves slower in the vicinity of the 
walls than in the center of the channels). This stretch- 
ing effect, is balanced by transverse molecular diffusion, 
resulting a Taylor- like dispersion mechanism (section [Tl|. 
This effect is made significant by the increase with Pe 
of the correlation length of dye transport along chains of 
flow channels parallel to the mean flow discussed above 
in sections HVEl and [IVll 

The effect of the local velocity gradients is also enhanced 
by the specific topology of 2D micromodels. The upper 
and lower walls are indeed continuous and some flow lines 
remain close to them over their full length: as in Taylor 
dispersion, slow solute particles near these walls may only 



move away from them through molecular diffusion. Simi- 
larly, fast moving particles half way between the walls can 
only reach them through transverse molecular diffusion. 
The large correlation length of these velocity contrasts 
also results in Taylor-like dispersion effects. 
Yet, the influence of the disorder of the medium cannot 
be completely neglected (since some tracer always moves 
into lower velocity transverse channels): the global dis- 
persion results therefore from the combined effects of ge- 
ometrical and Taylor dispersion. As a result, the macro- 
scopic dispersion coefficient D follows in this regime a 
power law D oc Pe^ of the Peclet number with an ex- 
ponent /3 ~ 1.3 intermediate between the values 1 and 
2 corresponding respectively to Taylor and geometrical 
dispersion. 

Regarding the influence of the shear thinning properties, 
increasing the polymer concentration enhances velocity 
contrasts between different flow channels while it flat- 
tens the velocity profile in individual channels. The first 
effect increases geometrical dispersion and is indeed ob- 
served at low Peclet numbers. The second reduces Taylor 
dispersion : this explains why, at higher Pe values, the 
values of Id for the two polymer solutions are similar when 
the influence of Taylor dispersion is large. 
To conclude, the dispersion measurements reported in 
the present work for transparent micromodels provide 
significant novel informations on the influence of the flow 
velocity and fluid rheology on miscible displacements in 
porous media. Quantitative high resolution optical mea- 
surements have allowed for thorough studies over a broad 
range of length scales: it has in particular been possible 
in the same experiment both to determine macroscopic 
parameters such as the effective dispersivity and to anal- 
yse at the pore scale the dynamics of concentration vari- 
ations. 

In particular, the local analysis of the time lag between 
the invasions of longitudinal and transverse channels of 
the model has allowed us to relate the transition between 
two dispersion regimes for Pe ~ 10 to variations of mix- 
ing in channel junctions. The variations of small scale 
structures of the displacement front with the Peclet num- 
ber and the polymer concentration also provides informa- 
tions on the spatial correlation of transport at the local 
scale. 

In the future, investigation of these effects at still higher 
resolutions should allow for detailed direct studies of mix- 
ing and flow patterns right inside individual flow channels 
and their junctions. 
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